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Abstract 

Miscanthus x giganteus is exceptional among C 4 plants in its ability to acclimate to chilling (<14 °C) and maintain a 
high photosynthetic capacity, in sharp contrast to maize, leading to very high productivity even in cool temperate cli- 
mates. To identify the mechanisms that underlie this acclimation, RNA was isolated from M x giganteus leaves in chill- 
ing and nonchilling conditions and hybridized to microarrays developed for its close relative Zea mays. Among 21 000 
array probes that yielded robust signals, 723 showed significant expression change under chilling. Approximately half 
of these were for annotated genes. Thirty genes associated with chloroplast membrane function were all upregulated. 
Increases in transcripts for the Ihcb5 (chlorophyll a/b-binding protein CP26), ndhF (NADH dehydrogenase F, chloro- 
plast), atpA (ATP synthase alpha subunit), psbA (D1), petA (cytochrome f), and Ihcb4 (chlorophyll a/b-binding protein 
CP29), relative to housekeeping genes in M. x giganteus, were confirmed by quantitative reverse-transcription PCR. 
In contrast, psbol, Ihcb5, psbA, and Ihcb4 were all significantly decreased in Z. mays after 14 days of chilling. Western 
blot analysis of the D1 protein and LHCII type II chlorophyll a/b-binding protein also showed significant increases in M. 
x giganteus during chilling and significant decreases in Z. mays. Compared to other C 4 species, M. x giganteus grown 
in chilling conditions appears to counteract the loss of photosynthetic proteins and proteins protecting photosystem 
II typically observed in other species by increasing mRNA levels for their synthesis. 

Keywords: C 4 photosynthesis, chilling, chlorophyll a/b-binding protein, cold, D1 protein, LHCII, low temperature, maize, 
Miscanthus, transcription. 



Introduction 

Miscanthus x giganteus (Greef & Deuter ex Hodkinson & 
Renvoize; Greef and Deuter, 1993; Hodkinson and Renvoize, 
2001), a rhizomatous perennial grass and bioenergy crop, 
appears unique among C 4 species in its productivity and 
maintenance of photosynthetic capacity during chilling 
(<14 °C) (Beale and Long, 1995; Beale et al, 1996; Naidu and 



Long, 2004; Long and Spence, 2013). In the closely related C 4 
crop Zea mays L., photosynthetic capacity declines rapidly 
in response to chilling (Long et al, 1983; Naidu et al, 2003; 
Wang et al, 2008a). This limits the growth of Z. mays in the 
US corn belt to the chilling-free period of the year, despite 
intensive breeding efforts to improve chilling tolerance 



© The Author 2014. Published by Oxford University Press on behalf of the Society for Experimental Biology. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.Org/licenses/by/3.0/), which 
permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. 



3738 | Spence ef a/. 



(Rodriguez et al, 2010). In a side-by-side comparison, M. x 
giganteus proved 59% more productive than Z. mays within 
the US corn belt even though both crops converted inter- 
cepted solar radiation into biomass with equal efficiency 
(Dohleman and Long, 2009). The higher productivity of M. 
x giganteus results from its ability to produce photosyntheti- 
cally active leaves and maintain them earlier and to retain 
photosynthetically competent leaves later in the growing sea- 
son — times when chilling temperatures damage photosynthe- 
sis in Z. mays (Dohleman and Long, 2009; Long and Spence, 
2013). 

Chilling during daylight reduces the maximum quan- 
tum yield of photosynthetic C0 2 assimilation (<£C0 2max ) hi 
Z. mays due to a combination of chilling-dependent photoin- 
hibition and impaired synthesis of key proteins of the PSII 
and the light-harvesting complex (LHC) (Nie and Baker, 
1991; Long et al, 1994; Fryer et al, 1995). M. x giganteus 
shows little depression in <5C(9 2max and photosynthetic rate 
during chilling (Beale et al., 1996; Naidu and Long, 2004; 
Farage et al, 2006; Wang et al, 2008a; Purdy et al, 2013). 
When grown side-by-side in controlled environments at 14 
°C, < t > C0 2m ^ was decreased significantly, relative to 25 °C, by 
c.50% in Z. mays, but was not significantly decreased in M. x 
giganteus. In the same study, photosystem II (PSII) maximum 
dark-adapted efficiency (F v /F m ) and operating efficiency 
( Opsn) hi saturating light were reduced c.20 and 80%, respec- 
tively, in Z. mays yet were hardly affected in M. x giganteus 
(Naidu and Long, 1994 J. Previous work has also compared 
the photosynthetic capacity of M. x giganteus under chill- 
ing conditions to a population of the C 4 sedge Cyperus lon- 
gus native to southern England at 52° N — an exceptionally 
high latitude for any C 4 species (Farage et al. , 2006; Long and 
Spence, 2013). Growth of M. x giganteus at chilling tempera- 
tures >8 °C did not alter the number of carbon dioxide mol- 
ecules fixed per electron transferred through PSII (A/J PSII ). In 
contrast, A/J PSII declined significantly in C. longus, suggest- 
ing that even among C 4 species native to cool climates, M x 
giganteus is exceptional in its ability to maintain high photo- 
synthetic efficiency (Farage et al, 2006). This maintenance 
of photosynthetic efficiency in M. x giganteus may in part 
result from avoidance of chilling-dependent photoinhibi- 
tion through the dissipation of excess absorbed light energy, 
termed nonphotochemical quenching (NPQ). Zeaxanthin, 
a photoprotective xanthophyll involved in inducible NPQ, 
increased 20-fold in M. x giganteus grown at 8 °C (Farage 
et al. , 2006). The maintenance of photosynthetic C0 2 assimi- 
lation by M. x giganteus would further dissipate absorbed 
excitation energy. Both factors are pivotal in the avoidance of 
chilling-dependent photoinhibition, which in Z. mays is char- 
acterized by the accumulation of nonfunctional PSII reaction 
centres (Fryer et al, 1995). 

M. x giganteus grown at 14 °C reaches nearly the same 
light-saturated rates of C0 2 assimilation (^4 sat ) as plants 
grown at 25 °C, when both are measured at 25 °C ('Naidu and 
Long., 2004; Wang et al, 2008a). Conversely, Z. mays leaves 
from plants grown at 14 °C display a light-saturated leaf C0 2 
uptake rate (A sat ) <10% of that of plants grown at 25 °C, 
when both are measured at 25 °C (Naidu and Long, 2004). 



This ability of M. x giganteus to maintain A slit under pro- 
longed chilling conditions has been correlated with its ability 
to maintain levels and activities of pyruvate phosphate diki- 
nase (PPDK) and Rubisco when grown at 14 °C that equal or 
exceed levels for plants at 25 °C (Wang et al, 2008a,b). These 
two enzymes have been shown to have high metabolic control 
coefficients for light-saturated C 4 photosynthesis (Furbank 
et al, 1997). In contrast, the concentrations and activities of 
these two enzymes decline markedly in Z. mays when trans- 
ferred to 14 °C (Naidu et al. , 2003; Wang et al. , 2008a). While 
these differences likely account for maintenance of A Silt during 
chilling acclimation, they do not explain the parallel mainte- 
nance of a high maximum quantum yield of leaf cPC0 2max 
(Naidu and Long, 1994). The initial slope of the response 
of leaf C0 2 assimilation (A) to photon flux (i.e. < t>C0 2m:iX ) 
provides an important measure of capacity for light-limited 
photosynthesis. In Z. mays, a decline in <£C0 2max is attributed 
to loss of efficiency of PSII, which is associated with impair- 
ment in the synthesis of key PSII and LHC proteins (Nie and 
Baker, 1991; Fryer et al, 1995; Caffarri et al, 2005/ 

While the comparative physiological aspects of chilling tol- 
erance of photosynthesis in M. x giganteus relative to Z. mays 
have been studied extensively, the molecular basis of this dif- 
ference is largely unexplored. Analysis of gene expression 
in chilling-intolerant Z. mays during acclimation to 14 °C 
showed significant downregulation in transcripts encod- 
ing chlorophyll a/b-binding proteins, Rubisco, and PPDK 
(Trzcinska-Danielewicz et al, 2009; Wang et al, 2008a,fe). 
The current study asked whether the demonstrated physi- 
ological tolerance of M. x giganteus to chilling compared 
to Z. mays corresponds to maintenance or upregulation of 
a wider range of genes encoding key aspects of the photo- 
synthetic apparatus, using the chilling treatments of previous 
physiological studies (Wang et al, 2008a,o). Particular atten- 
tion is given to those proteins shown to become damaged or 
deficient in Z. mays during acclimation to chilling associated 
with the decline in <5C0 2max . 

Materials and methods 

Biological materials 

Plant propagation, chilling treatment, and leaf sample collection 
were performed in two independent experiments for M. x giganteus. 
The leaf samples from the first experiment were used in the micro- 
array analysis and quantitative reverse-transcription PCR (qPCR) 
validation, and the leaf samples from the second experiment, which 
included a side-by-side comparison with Z. mays, were used for 
the comparisons of transcripts by qPCR and amounts of specific 
proteins. 

Experiments followed the treatment protocol of Wang et al. 
(2008fl,/j). Briefly, plants were first grown in warm conditions of 
25 °C days and then transferred to 14 °C for 14 days. Prior studies 
have shown that this is sufficient time for completion of acclima- 
tion of photosynthesis in mature leaves of both species. In the first 
experiment, M. x giganteus rhizomes were propagated in soil-less 
potting media (Sunshine Mix LCI 2 , Sun Gro Horticulture, Bellevue, 
WA, USA) in 1 .2-1 pots following the procedure of Naidu and Long 
(2004). Eight individual plants were grown in two identical con- 
trolled-environment growth chambers (Conviron El 5, Controlled 
Environments, Winnipeg, Manitoba, Canada) under 500 umol 
photons m~ 2 s -1 incident light provided by an equal mixture of 
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high-pressure mercury and sodium lamps, in a 14/10 light/dark cycle 
at either 14/12 °C (chilling) or 25/20 °C (control). Both chambers 
were maintained at 70% relative humidity. The eight plants were ini- 
tially grown in the control chamber for 2 weeks, and then half of 
these plants, selected by fully randomized design, were transferred 
to the chilling conditions. After 2 weeks of chilling, and 1 month of 
total growth, leaf samples were taken 2-3 h after illumination began 
(i.e. mid-morning) from the midpoint of the youngest fully expanded 
leaf (Naidu and Long, 2004; Wang et al, 2008/?). Full expansion 
was defined by emergence of the ligule. Shoots had approximately 
four leaves each at the time of sample collection. All leaf samples 
were immediately plunged into liquid nitrogen and then transferred 
to storage at -80 °C until RNA extraction for microarray and qPCR 
analysis. 

The conditions for the second experiment were identical to the 
first apart from the following addition. Z. mays cv. B73 seeds were 
germinated in 1 .2 liter pots in soil-less potting media following the 
procedure of Naidu and Long (2004) and grown alongside M. x 
giganteus following the sampling and conditions described above. 



4 x 44K Agilent oligonucleotide two-colour cDNA maize 
microarray 

RNA extraction Total RNA was extracted from the leaf samples. 
Samples were homogenized in 1 ml TRIZOL Reagent (Invitrogen, 
Carlsbad, CA, USA) using a glass-Teflon homogenizer. The homog- 
enized samples were then incubated for 5 min at room temperature 
to permit the complete dissociation of the nucleoprotein complexes. 
Each 1ml sample had 0.2 ml chloroform added, was incubated at 
room temperature for 150 s, and was centrifuged at 12 000 g for 
1 5 min at 4 °C. The RNA in the resulting supernatant was precip- 
itated with 0.5 ml isopropyl alcohol and was isolated by centrifu- 
gation at 12 000 g for 10 min at 4 °C. The RNA pellet was dried 
and redissolved in 50-jj.l RNase-free water for 10 min at 55 °C and 
then further purified with the RNeasy Kit (Qiagen, USA). Final 
RNA quantity was determined spectrophotometrically at X 260 nm 
(NanoDrop ND1000, NanoDrop Technologies, Wilmington, DE, 
USA). The purity of RNA was assessed by the 260/280 absorbance 
ratio. Only those RNA samples with ratios of >2.0, indicating no 
significant presence of contaminants, were used. 
Microarray fluorescent cDNA synthesis and hybridization cDNA 
was synthesized from total RNA using a Agilent Low RNA Input 
Linear Amplification Kit, with spike-in controls (5184-3523), fol- 
lowing the manufacturer's instructions (Agilent Technologies, Santa 
Clara, CA, USA). Two of the four control samples and two of the 
four treatment samples were labelled with Cy3 dye and the other 
samples were labelled with Cy5 dye, to account for any possible dye- 
bias that may occur during hybridization or fluorescence detection 
during scanning of the slides. Agilent 4 x 44K oligonucleotide two- 
colour cDNA Maize microarray slides were used. Hybridization 
and washing was performed following the Agilent Two-Colour 
Microarray-Based Gene Expression Analysis protocol (Agilent 
Technologies). A total of four 4 x 44k slides (i.e. four technical repli- 
cates) were used for each of the four biological replicates. Hybridized 
slides were scanned with a GenePix4000B microarray scanner (Axon 
Instruments, Concord, ON, Canada) and features were extracted 
and analysed using GenePix Pro 6. 1 (Axon Instruments). Each array 
was inspected individually and low-quality spots (i.e. those with flu- 
orescence levels that were not significantly distinguishable from the 
background fluorescence in the red-to-green channel intensity ratio 
histograms of the array) were manually flagged and eliminated from 
further analysis. 

Microarray statistical analysis All calculations were performed 
using the R computing environment (Venables et al, 2013). 
Normalization and statistical analysis was done using the R package 
Limma provided by the Bioconductor repository (http://bioconduc- 
tor.org; Wettenhall and Smyth, 2004). No background correction 
was performed, a common practice for spotted arrays with low 
background (Smyth et al, 2003, 2005). The global loess function 



was employed for within-array normalization and the scale method 
was applied for between-array normalization (Smyth and Speed, 
2003). Test-statistics were determined using eBayes fit and topTable 
functions (Smyth, 2004). Use of the duplicate correlation function 
accounted for duplicated probe spots on the Agilent array and for 
technical replication (Smyth et al, 2005). The lmFit function was 
applied to fit the data to a linear model using quantitative nonnega- 
tive weights for flagged spots, and a design matrix was used to com- 
pare expression values in the control vs. treatment arrays. Genes that 
were differentially expressed at 14 °C versus 25 °C were defined on 
the basis of a significant P- value, adjusted for multiple testing, using 
FDR (adj. P-value 0.00001) and a log 2 fold-change cut off of 0.70 
or -0.70. These criteria were validated by the ability to also detect 
significant expression differences by qPCR, as defined previously 
(Wurmbach et al, 2003). 

Gene annotation and MapMan application 

Sequences for the FASTA 60-mer probes on the Agilent 4x44K 
oligonucleotide two-colour cDNA Maize microarray slides were 
provided by Dr. Virginia Walbot at Stanford University. A gene 
annotation file for the Agilent 60-mer probe sequences was gener- 
ated by BLASTN of probe sequences to the Z. mays B73 reference 
genome sequence, build 2 (Schnable et al, 2009). For genes found to 
change significantly between control and treatment, but lacking any 
published annotation at the time of the experiment, annotation was 
attempted by performing a tblastn BLAST query (http://www.ncbi. 
nlm.nih.gov/) on the corresponding 60-mer probe FASTA sequence. 

To visualize the differentially expressed genes, the MapMan soft- 
ware tool was used (Usadel et al, 2009). The MapMan software uti- 
lizes the Affymetrix platform to generate hierarchical categories of 
gene function and pathway placement. In order to use this software, 
to visualize results obtained here on the Agilent microarray, a 'map- 
ping file' was created based upon the pre-existing Z. mays ontology/ 
mapping file with Affymetrix-assigned probe identifiers. The Agilent 
annotations were cross-referenced with the Z. mays mapping file 
and gene annotations matching those in the Z. mays mapping file 
BIN or sub-BIN for functional categories. The Affymetrix identifier 
was replaced with the Agilent identifier. The Agilent identifier cor- 
responding to the assigned BIN/category was then used to create an 
experiment file containing the corresponding down- or upregulated 
log 2 fold-change value. This experiment file and the new mapping 
file was then uploaded into the MapMan Image Annotator module 
(Usadel et al, 2009). 



Quantitative reverse-transcription PCR 

Quantitative reverse-transcription PCR was used to validate the 
up- and downregulation of transcripts of particular interest that 
were found to be differentially expressed in the microarray experi- 
ment. Genes for validation with known Z. mays sequences (http:// 
www.ncbi.nlm.nih.gov/) were compared against available M. x 
giganteus genome survey sequences obtained by 454 sequencing 
(Swaminathan et al, 2010) and transcripts assembled from Illumina 
short-read sequencing (Barling et al, 2013). Sequences with more 
than 98% identity were candidates for qPCR validation. Six tran- 
scripts were chosen: ndhF (NADH dehydrogenase F, chloroplast), 
atpA (ATP synthase alpha subunit), IhcbS (chlorophyll a/b-binding 
protein CP26), psbol (oxygen-evolving enhancer protein 1), aps 
(ATP synthase), and tps (terpene synthase). The primers for the tran- 
scripts were designed based on M. x giganteus sequences (Table 1). 
Four candidate reference genes for qPCR were chosen because they 
are common housekeeping genes that have stable expression levels 
between the control and treatment samples in the microarray experi- 
ment: tuba6 (alpha-tubulin 6), tubb3 and tubb6 (beta-tubulin 3 and 
6), and tubgl (gamma-tubulin 1). 

SuperScript-III RT (Invitrogen) was used to synthesize cDNA 
from the same M. x giganteus samples used for the microarray 
experiment. SYBR Green Master Mix (Applied Biosystems, Foster 
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Table 1 . Primer sequences for qPCR validation of microarray results and comparison with Z. mays 



Transcript 



Forward (5'-3') 



Reverse (5 -3 ) 



Expected product size (bp) 



tuba6 

tubb3 

tubbS 

tubb6 

tubgl 

act 

ndhF 

atpA 

Ihcb4 

Ihcb5 

psbol 

aps 

tps 

psbA 

petA 



AAGACGCAGCTAACAACTTTGCCC 

ATGGACGAGATGGAGTTCACCGAA 

CTTTGTTTCACCTGCACCGCTTGA 

TTCTGACCTTCAGTTGGAGCGTGT 

ATGCACCTGAGTGGCGAAGTCCT 

TGAGGCCACGTACAACTCCATCAT 

ACCCACTCCCATTTCGGCTCTTAT 

ATGGCGGATTCACCCGCTACATTA 

AGTAGAAGATGAGCATGGCGAGCA 

TCGCCATGTTCTCCATGTTCGGAT 

ACCTACGACGAGATCCAGAGCAA 

TGGATGGAGCAGAAGAACAGGGTT 

ATGCGAGCATCAGACAGTTCAGGA 

ACTTAGTTTCCGTCTGGGTATGCG 

CGCACATCTATTTCAAACGCATA 



CACCACCAACAGCATTGAACACCA 

TGATGTACCACAGCAGCGAAGACT 

AGGAAGGAACCATCACAGGAGCAA 

TGCCCAAACACAAAGTTGTCAGGG 

ATCGGGATGGAGTTCTGGAAGCA 

CCTTTCAGGTGGCGCAATCACTTT 

GAGCAAGAGCTAAAGTGGCTCCTA 

TGGCGATAAGCTTGTGCCTGTTTG 

AGAACTTCGCCAACTTCACCG 

TAGGCCCAGGCGTTGTTGTTGA 

GCACAGCTTCTTCAGCTGGTACTT 

TGCGGTATTTGCCTTTCAGCTTCG 

CAGCTACCATTTGCACGGCAGAAA 

TAAGGATGTTGTGCTCTGCCTGGA 

TACAATTCGTCCAGTTGCTTCCCG 



140 
174 

98 
164 

81 
180 
166 
136 
180 
114 
134 
120 

88 
185 

81 



City, CA, USA) was used for qPCR and samples were amplified on a 
high-throughput microwell-plate-based thermo-cycler (LightCycler 
480 System, Hoffmann-La Roche, Indianapolis, IN, USA). Results 
were normalized against tubb6 and analysed using the AACT method 
(Schmittgen and Livak, 2008). tubb6 was chosen for standardization 
as it exhibited the lowest variance across all samples included in the 
qPCR experiment. One-tailed Student's t-test was used to determine 
a significant increase in transcript level at a<0.05 due to chilling for 
each gene. 

The second experiment comparing M. x giganteus to Z. mays 
grown side by side in control and chilling conditions used the same 
protocol. Three additional transcripts were added with known 
involvement in chilling response in Z. mays: psbA (Dl), pet A 
(cytochrome f), and lhcb4 (chlorophyll a/b-binding protein CP29) 
(Trzcinska-Danielewicz et ah, 2009; Table 1). Act was used as the 
reference gene for normalization of qPCR results because it shows 
stable expression in both species and has been previously shown 
to have stable expression during chilling in Z. mays (Trzcinska- 
Danielewicz et al. , 2009). Results were normalized against Act and 
analysed using the AACT method (Schmittgen and Livak, 2008). 
Student's t-test was used to determine significant differences at 
a<0.05 between treatments for each species and gene. 



Western blot analysis 

Total leaf protein was extracted from leaf discs of area 1.27 cm 2 
punched from the centre of the lamina of each biological replicate 
(n=4) of each species. Proteins were extracted by grinding frozen 
samples in individual tubes with stainless steel grinding balls (GBSS 
156-5000-01, OPS Diagnostics, Lebanon, NJ, USA) twice for 30 
s at 300 strokes min 1 (2000 Geno/Grinder, SPEX SamplePrep, 
Metuchen, NJ, USA). The samples were placed on dry ice and 500 ul 
extraction buffer added to each tube, as described by Wang et al. 
(20086). Briefly, the extraction buffer contained 50 mM HEPES 
(pH 8.0) 0.05% (v/v) Triton-X 100, ImM EDTA (pH 8.0), lOmM 
MgCl 2 , 5.96 mM dithiothreitol, 1% (w/v) casein, 1% (w/v) polyvi- 
nylpyrrolidone, and one protease inhibitor cocktail tablet per 10 ml 
of extraction buffer (cOmplete ULTRA Tablets, Mini, EDTA-free 
tablets, Roche Applied Science, Indianapolis, IN, USA). The ground 
tissue and buffer was homogenized and centrifuged for lmin at 15 
000 g and then the resulting supernatant was stored at -80 °C until 
Western blot analysis. 

Total leaf protein, on an equal leaf area basis, was loaded onto 
a 10-lane SDS-PAGE 10% Tris-HCL gel, such that each lane rep- 
resented the same area of leaf (Mini-PROTEAN TGX 456-1033, 
Bio-Rad Life Science, Hercules, CA, USA). A leaf area basis, rather 
than the usual total protein basis, was used so that findings could 



be related to changes in photosynthetic capacity, which is measured 
on a leaf area basis (Naidu and Long, 1994; Wang et al, 20086). 
Polypeptides were separated by SDS-PAGE, as described previously 
(Wang et al. , 20086). Each gel contained four biological replicates 
for the control and of the chilling treatment for both M. x gigan- 
teus and Z. mays. Protein-prestained standards were used to pro- 
vide a molecular weight ladder (Precision Plus Protein Kaleidoscope 
Standards 161-0375, Bio-Rad Life Science). Separated proteins 
were blotted onto polyvinylidene fluoride membranes in transfer 
buffer (0.025 M Tris base, 0.192M glycine, and 20%, v/v, methanol) 
at 100 V for lh at 4 °C. Membranes were then incubated with pri- 
mary polyclonal antibodies in Tris-buffered saline (170-6435, Bio- 
Rad Life Science) with 0.05% Tween 20 (P5927, Sigma-Aldrich, St. 
Louis, MO, USA). LHCII chlorophyll a/b-binding proteins and 
the Dl protein and its degradation products were detected with 
rabbit polyclonal antibodies (AS01-003 and AS05-084, Agrisera 
Antibodies, Vannas, Sweden). Membranes were then incubated with 
anti-rabbit IgG (whole molecule) alkaline phosphatase antibody 
produced in goat (A7539, Sigma-Aldrich). The secondary antibody 
was detected with a mixture of 5-bromo-4-chloro-3-indolyl-phos- 
phate and nitroblue tetrazolium, which binds alkaline phosphatase 
(Western Blue Stabilized Substrate, Promega, Madison, WI, USA). 
Digital images were made of each membrane and the relative vol- 
ume and intensity for each band was quantified using the 2D image 
analysis software ImageQuant TL 7.0 (GE Healthcare Biosciences, 
Little Chalfont, UK). 

Mixed-model analysis of variance followed by least significant 
difference tests of differences between means determined the sig- 
nificance of the effect of chilling on the content of each of LHCII 
chlorophyll a/b-binding protein, Dl protein, and Dl degradation 
products in M. x giganteus and Z. mays (PROC MIXED SAS ver- 
sion 9.2, SAS Institute, Cary, NC, USA). Results were considered 
significant at a<0.05 (n-4). 



Results 

M. x giganteus leaf cDNA hybridized to the Agilent 44K 
maize oligonucelotide microarray with 48% efficiency. 
Approximately 21 000 probes showed a median fluorescence 
intensity that could be clearly distinguished from the back- 
ground. Chilling changed the relative signal intensity for 723 
probes relative to the control; based on a significance thresh- 
old of a log 2 fold-change of >0.7 and false-discovery-rate- 
adjusted probability of a<0.00001. These thresholds were 
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selected to include candidate genes with small but reproduc- 
ible changes in expression (Wurmbach et al, 2003). Chilling 
caused significant increase in signal intensity for 410 probes 
and a decrease in 3 1 3 . Approximately half of the probes show- 
ing significant changes corresponded to sequences annotated 
for predicted or putative functions in Z. mays, and these were 
categorized accordingly in MapMan (Fig. 1). A complete list 
of putative chilling-responsive genes identified here is given in 
the supporting information (Supplementary Tables SI and S2 
available at JXB online). Functional categories that were 
significantly represented among chilling-responsive genes 
included photosynthesis, stress response, RNA processing, 
protein synthesis, targeting and degradation, transport, sec- 
ondary metabolism, and hormone metabolism. 

Most striking among these categories were genes associ- 
ated with the light reactions. Among the probes binned by 
MapMan to this category, 30 were significantly upregu- 
lated and none were significantly downregulated (Table 2). 
Particularly prominent were transcripts for PSII-associated 
proteins (Fig. 2 and Table 2). Within the carbon metabolism 
pathway, the probes for PPDK and Rubisco were also main- 
tained and showed high expression values (Supplementary 
Table S3). Probes for real and real, the genes for Rubisco 
activase, an enzyme that is critical to the activation and sta- 
bility of Rubisco (Portis, 2003), were upregulated 145%, 
although just below the fold-change cut off (log 2 < 0.7). 

Six of the probes found to change on the microarray were 
selected for validation by qPCR in the first experiment. 
Each showed the same directional change as well as similar 
quantitative changes in expression to that observed from 
the microarray (a<0.05; Fig. 3). These included four upreg- 
ulated probes representing genes encoding key proteins 



of photosynthetic light reactions: psbol (oxygen-evolving 
enhancer protein 1), lhcb5 (chlorophyll a/b-binding protein 
CP26), ndhF (chloroplast NADH dehydrogenase subunit 
5), and atpA (ATP synthase alpha subunit), as well as two 
separate downregulated probes, one for aps (ATP sulphury- 
lase) and the other for tps (terpene synthase). Five probes 
representing genes with likely housekeeping functions were 
also tested for constancy of expression between the treat- 
ment and control. tubb6 had the most stable expression 
among these genes and was therefore used as a reference 
gene for each of the other seven transcripts in estimating 
relative fold-changes by the AACT method (Schmittgen 
and Livak, 2008). 

In a second qPCR experiment comparing M. x giganteus 
to Z. mays, using five transcripts validated from the microar- 
ray experiment and three additional genes with known regula- 
tion during chilling, significant differences in gene expression 
between the two species were found during chilling. Changes 
in transcripts for M. x giganteus were consistent with those 
found in the first experiment, although psbol expression 
change was not significantly different (P=0.14). In Z. mays, 
transcripts for psbA (Dl), lhcb4 (chlorophyll a/b-binding pro- 
tein CP29), psbol, and lhcb5 all showed significant downregu- 
lation (a<0.05), while the same transcripts, with the exception 
of psbol, in M. x giganteus were significantly upregulated 
(a<0.05) (Fig. A), pet A (cytochrome F) indicated downregula- 
tion in Z. mays but it was not significantly different from con- 
trol (P=0.14). Similar expression patterns were seen in ndhF 
for both species (Fig. 4). These two qPCR experiments were 
performed on RNA samples from completely independent 
chilling treatments and showed the apparent reproducibility 
of these results in M. x giganteus. 
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Fig. 1. Functional categorization of transcripts found to be up- or downregulated during chilling (<14 °C) in Miscanthus x giganteus. Functional groups 
are as categorized by MapMan (Usadel era/. 2009). Red boxes indicate downregulated transcripts and blue boxes represent upregulated transcripts, as 
indicated by the scale from -1 .6 to 1 .6 log 2 fold-change in expression. See Table 2 for numerical details of the change in each transcript. 
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Table 2. Relative changes in abundance of all M. x giganteus leaf transcripts for chloroplast proteins, associated with the light reactions 
that showed a log 2 fold increase >0. 70 with chilling (14 °C) 

Data are relative to control (25 °C). There were no significant decreases in response to chilling. 



Transcript 


Functional category 


I0Q2 fo 


Chlorophyll a/b-binding apoprotein CP26 precursor 


Photosystem II 


light harvesting complex II 


0.815 


Photosystem II P680 chlorophyll A apoprotein (CP-47 protein) 


Photosystem II 


light harvesting complex II 


0.93 


Violaxanthin de-epoxidase precursor 


Photosystem II 


light harvesting complex II 


0.855 


D1 protease precursor (fragment) 


Photosystem II 


light harvesting complex II 


1.751 


Photosystem II 44kDa reaction centre protein (P6 protein) (CP43) 


Photosystem II 


reaction centre proteins 


0.95 


Photosystem II reaction centre J protein 


Photosystem II 


reaction centre proteins 


0.735 


Oxygen evolving enhancer protein 3 


Photosystem II 


oxygen-evolving proteins 


0.989 


Thylakoid lumenal 17.4kDa protein, chloroplast 


Photosystem II 


thylakoid lumenal proteins 


1.408 


Thylakoid lumenal 21 .5kDa protein, chloroplast precursor 


Photosystem II 


thylakoid lumenal proteins 


0.805 


Thylakoid lumenal 15kDa protein, chloroplast precursor (p15) 


Photosystem II 


thylakoid lumenal proteins 


0.837 


NDA2 (ALTERNATIVE NAD(P)H DEHYDROGENASE 2); NADH 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


0.906 


dehydrogenase 








NAD(P)H-quinoneoxidoreductase chain 5, chloroplast (NADH- 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


1.081 


plastoquinoneoxidoreductase chain 5) 








NAD(P)H-quinoneoxidoreductase chain 1 , chloroplast (NADH- 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


1.914 


plastoquinoneoxidoreductase chain 1) 








NAD(P)H-quinoneoxidoreductase chain 4L, chloroplast (NADH- 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


1.052 


plastoquinoneoxidoreductase chain 4L) 








NADH-plastoquinoneoxidoreductase subunit 5 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


0.856 


NADH-plastoquinoneoxidoreductase subunit 4 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


0.766 


NAD(P)H-quinoneoxidoreductase chain 1 , chloroplast (NADH- 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


0.797 


plastoquinoneoxidoreductase chain 1) 








Cytochrome b6 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


1.223 


3Fe-4S ferredoxin 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


2.734 


Ferredoxin-nitrite reductase precursor-maize (fragment) (Zea mays) 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


1.153 


Ferredoxin-NADP reductase, root isozyme, chloroplast precursor (FNR) 


Photosystem II 


NAD(P)H dehydrogenase/PQ 


1.207 


ATP synthase alpha chain 


Photosystem light reactions: ATP synthase 


1.007 


Cytochrome c biogenesis protein family 


Chloroplast biogenesis 


0.972 


Chloroplast ORF70 


Chloroplast biogenesis 


0.786 


matK protein (trnK intron) 


Chloroplast proteins 


0.8146 


Zea mays chloroplast rRNA-operon 


Chloroplast RNA regulation of transcription 


1.078 


50S ribosomal protein L29, chloroplast precursor 


Chloroplast protein synthesis: plastid ribosomal proteins 


0.907 


Chloroplast 30S ribosomal protein S18 


Chloroplast protein synthesis: plastid ribosomal proteins 


0.972 


50S ribosomal protein L21 , chloroplast precursor (CL21) 


Chloroplast protein synthesis: plastid ribosomal proteins 


0.721 


Chloroplast 50S ribosomal protein L1 6 


Chloroplast protein synthesis: plastid ribosomal proteins 


0.775 



The second experiment also showed that the pattern of 
change in transcripts was reflected in the amounts of two key 
proteins. As predicted from the transcript changes, amounts 
of LHCII type II chlorophyll a/b-binding protein and Dl 
protein per unit leaf area increased significantly (a<0.05) in 
M. x giganteus leaves in response to chilling, but decreased 
significantly (a<0.05) in Z. mays (Fig. 5). The proportion 
of Dl degradation products relative to Dl protein did not 
change significantly in either species (Fig. 5). 

Discussion 

The present study asked whether the physiological chilling 
tolerance of photosynthesis in M. x giganteus corresponded 
to maintenance or upregulation of genes encoding key 
aspects of the photosynthetic apparatus, specifically those 
proteins shown to be damaged or deficient in Z. mays dur- 
ing acclimation to chilling that relate to energy transduction 
on the photosynthetic membrane. Decline in the maximum 



quantum yield of <PC0 2max in maize with chilling has been 
associated with the loss of key proteins of the photosynthetic 
membrane, particularly those of PSII and those involved 
in transfer of excitation energy to PSII (Fryer et al, 1995). 
In contrast, a small but not statistically significant loss of 
<r , C0 2max was indicated when M. x giganteus was subjected 
to the same chilling conditions (Naidu and Long, 2004). This 
raises the question of how this species avoids this loss. Here, 
30 genes annotated as encoding proteins with key functions 
in the light reactions of photosynthesis showed significantly 
increased expression, and none of this category showed a 
decrease. This suggests that M. x giganteus may avoid the 
losses observed in maize, at least in part, by increasing synthe- 
sis of these key proteins. Across all of the MapMan categories, 
photosynthesis was the only one that did not show any down- 
regulated probes (Tables 2 and 3 and Fig. 2). Functions of 
the proteins encoded by upregulated transcripts include elec- 
tron transport, light harvesting, photosystem reaction cen- 
tres, ATP synthesis, and the xanthophyll cycle which protects 
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Fig. 2. Transcripts related to photosynthesis showing change in A/7, x giganteus after 14 days of chilling. Colouration scale as in Fig. 1 . See Table 2 for 
numerical details of the change in each transcript. 



PSII against photoinhibitory damage (Demmig-Adams and 
Adams, 1992; Fryer et al, 1995; Farage et al, 2006). In the 
current study, in Z. mays following 28 h of 14 °C chilling, four 
genes associated with the light reactions of photosynthesis 
were significantly downregulated. These encode the chloro- 
phyll a/b-binding protein precursor, chlorophyll a/b-binding 
apoprotein and minor antenna complex CP24 precursor, chlo- 
rophyll a/b-binding and minor antenna complex CP26, and 
NAD kinase chloroplast precursor (Trzcinska-Danielewicz 
et al, 2009). The current study, which used a longer chilling 
treatment, similarly showed that genes encoding key proteins 
of the photosynthetic light reactions (i.e. Dl, cytochrome f, 
PSII CP29, oxygen-evolving enhance protein 1, and chloro- 
phyll a/b-binding proteins) were significantly downregulated 
in maize but upregulated in M. x giganteus (Fig. 4). 

Parallel to these transcript studies, proteomic analysis of 
Z. mays leaves developed at 13 °C show large decreases in the 
minor antenna complexes CP24, CP26, and CP29 (Caffarri 
et al, 2005). In a similar experiment, Nie and Baker (1991) 
showed large decreases at 12 °C also in the 26-kDa LHCII 
apoprotein, as well as the 32-kDa Dl reaction centre protein, 
cytochrome f ( 3 3 kDa ) , cytochrome b6/f subunit I V ( 1 7 kDa) , 
and the a and |3 subunits of the coupling factor (58 and 
57 kDa). Consistent with these prior studies, Western blotting 
in the current study similarly showed significant decreases in 



both Dl reaction centre protein and LHCII chlorophyll a/b- 
binding protein in Z. mays after 14 days of chilling. In con- 
trast, these same proteins increased in M. x giganteus (Fig. 5). 
For plants that acclimate to chilling conditions, increased 
overall leaf protein contents are commonly observed. This 
study group has previously shown that acclimation to the 
chilling conditions used here lowers total protein per unit leaf 
area by 17% in Z. mays and increases content by 2% in M. x 
giganteus (Naidu et al, 2003). This suggests that the much 
larger increases observed here in LHCII chlorophyll a/b- 
binding protein and Dl are more than the result of an overall 
increase in protein content. This difference in levels of these 
light reaction-related transcripts and protein contents is likely 
to be a key underlying mechanism, explaining an unusual but 
effective acclimation of photosynthesis to chilling. 

A critical function necessary to the maintenance of pho- 
tosynthesis during chilling in C 3 and C 4 plants is the ability 
of the chloroplast to degrade damaged Dl, a component of 
PSII, and to then synthesize and assemble Dl back into PSII 
(Bredenkamp and Baker, 1994; Fryer et al, 1995; Ruelland 
et al, 2009). This repair function is inhibited in Z. mays dur- 
ing chilling and this has been attributed to decreased expres- 
sion of psbA, the plastid gene that encodes Dl (Bredenkamp 
and Baker, 1994; Allen and Ort, 2001). The current results 
show that the transcripts for the gene for Dl protease, 
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Fig. 3. qPCR validation of six transcripts found to be up- or downregulated on the microarray (Table 2): psbOl (oxygen-evolving enhancer protein 1), 
Ihcb5 (chlorophyll a/b-binding protein CP26), ndhF (NADH dehydrogenase F, chloroplast, atpA (ATP synthase alpha subunit), aps (ATP sulphurylase), 
and tps (terpene synthase). The AACT method was used to calculate fold-change and all results were normalized against the reference gene tubb6. 
Grey bars indicate the qPCR-validated fold-change (primary Y-axis) and black squares indicate the microarray fold-change (secondary Y-axis). Error bars 
represent the standard error of the mean of the biological replicates (n=4). 
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Fig. 4. Comparison of gene expression between M. x giganteus and Z. mays after 14 days of chilling (14 °C) versus control (25 °C). qPCR-tested 
transcripts included five transcripts found to be up- or downregulated on the microarray (atpA, ATP synthase alpha subunit; psbol , oxygen-evolving 
enhancer protein 1 ; Ihcb5, chlorophyll a/b-binding protein CP26; ndhF, NADH dehydrogenase F; and tps, terpene synthase), and three additional light 
reaction transcripts: psbA, D1 ; petA, cytochrome f; and Ihcb4, chlorophyll a/b-binding protein CP29). Results were normalized against the reference 
gene act and the AACT method was used to calculate fold-change. Error bars represent the standard error of the mean of all biological replicates (n=4). 
Asterisks indicate significant differential expression between control and chilling (a< 0.05). 



involved in the degradation of damaged Dl , and for psbA, the 
gene encoding Dl, are strongly upregulated in M. x giganteus 
grown under chilling conditions, with fold-changes of 3.37 
(Table 2) and 2.06 (Fig. 4), respectively. Damage to Dl is a 
major cause of loss of PSII efficiency, since absorbed light 
energy will be channelled to inactive centres that, under light- 
limiting conditions, will lower <£C0 2 max (reviewed by Long 
et al, 1994). Previously a high expression of xanthophyll 



cycle pigments and accumulation of zeaxanthin during chill- 
ing of M. x giganteus has been shown (Farage et al, 2006). 
The present study provides an explanation for this: a signifi- 
cant increase in expression of the transcript for a violaxan- 
thin de-epoxidase, which converts violaxanthin to zeaxanthin 
(log 2 fold-change 0.855, Table 2). Accumulation of zeaxan- 
thin would serve to lessen damage to Dl by dissipating excess 
excitation energy by thermal de-excitation (Farage et al, 



Chilling tolerance in Miscanthus x giganteus | 3745 





M. x giganteus 


Z. mays 












25 kD 






LHCIItypall Chi a/b 


' 




binding protein 










32 kD 






D1 


25 kD 






D1 D»gr»d»d 











Species 


LHCII % 
Change 


a 


Dl % 
Change 


a 


Dl Degraded 
% Change 


a 


M. x g. 


+21.93 (±7.53) 


0.0311 


+20.69 (±6.52) 


0.0350 


-6.76 (±6.52) 


0.4845 


Z. mays 


-62.67 (±6.52) 


0.0001 


-31.26 (±6.52) 


0.0018 


-17.99 (±6.52) 


0.0657 



Fig. 5. Representative Western blot images showing changes in LHCII type II chlorophyll a/b-binding protein, D1 protein, and D1 degradation products 
after 1 4 days of chilling (1 4 °C) compared to control (25 °C) in M. x giganteus and Z. mays. The same leaves as for Fig. 4 were used. Values are 
mean±SE for all biological samples. Differences between control and chilling were assessed by Student's t-test (a<0.05, n=4). 

Table 3. Comparison of Miscanthus x giganteus and Zea mays transcriptional response to chilling 

1 indicates upregulation during chilling; J. indicates downregulation during chilling. Arrows in both directions indicate that some transcripts in this 
category were upregulated and others downregulated. A larger arrow indicates that more transcripts were up- or downregulated, and vice- 
versa. Transcriptional data: M. x giganteus, this study; Z. mays, Trzcinska-Danielewicz et al. (2009). 



Category 


M. x giganteus 




Z. mays 




Response 


No. of genes 


Response 


No. of genes 


Protein synthesis 


Ti 


8 


Ti 


8 


Protein degradation 


Ti 


17 


T 


1 


Chlorophyll biosynthesis 


T 


3 


I 


1 


PS II reaction centre and LHCII 


T 


6 


I 


3 


Temperature Induced Proteins 


I 


1 




0 


Cellulose synthesis 


I 


2 


I 


1 


Flavonol biosynthesis 


T4 


5 




0 


Brassinosteroid biosynthesis 


1 


1 




0 


Ethylene 


1 


1 


I 


3 


Gibberellic acid biosynthesis 


I 


3 




0 


Abscisic acid biosynthesis 


iT 


3 


I 


1 


Cytokinins 


1 


1 




0 


Wax biosynthesis 


TI 


2 




0 



2006). This, combined with increased capacity for processing 
damaged Dl should greatly decrease the number of centres 
with damaged Dl at any point in time and could explain the 
maintenance of high <J>C0 2max during chilling, as observed 
in M. x giganteus in chilling in contrast to maize (Naidu and 
Long, 1994; Beale et al, 1996). 

Maximum quantum yield ( < PC0 1 mnx) and light-saturated 
rates G4 8a t) of carbon assimilation are equally important in 
maintaining crop canopy carbon gain over a diurnal course, 
since this depends on both sunlit and shaded canopy leaves 
(Baker et al, 1988). While <PC0 2max in C 4 photosynthesis 
depends on the efficiency of energy transduction on the 
photosynthetic membrane and the stability of key proteins, 
such as Dl, A sllt depends on the maximum capacity of the 



rate-limiting steps in electron flow and carbon metabolism 
beyond the reaction centres, and diversion of electrons to 
sinks other than C0 2 assimilation. During chilling, M. x 
giganteus, in sharp contrast to Z. mays, maintains high lev- 
els of the enzymes of photosynthetic carbon metabolism 
that exert maximal metabolic control of A m , Rubisco, and 
PPDK. Rubisco protein contents are stable in M. x gigan- 
teus in response to chilling while they decline by more than 
50% in Z. mays during chilling (Naidu et al, 2003; Wang 
et al, 2008a). PPDK declines to an even greater extent in 
Z. mays but increases in M. x giganteus (Naidu et al, 2003; 
Wang et al, 2008/?). The current results correspond with 
these findings, with the expression levels of both Rubisco 
(log 2 fold-change 0.142) and PPDK (log 2 fold-change 0.015) 
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remaining stable during chilling. In addition, the transcripts 
for the two isoforms of Rubisco activase (real and real) 
were upregulated during chilling (Supplementary Table 
S3). Rubisco activase is a catalytic chaperone that is criti- 
cal to maintaining the stability and high catalytic activity of 
Rubisco. Specifically, Rubisco activase promotes the dissoci- 
ation of a wide variety of inhibitory sugar phosphates from 
the Rubisco active site (Portis, 2003). More recently, Shan 
et al. (2011) found that downregulation of rca was a key 
and necessary component in the jasmonate-induced senes- 
cence that is associated with a wide range of abiotic and 
biotic stresses. It is intriguing that chilling in M. x gigan- 
teus induced the opposite response to that seen in induced 
senescence. 

The Agilent 44K maize oligonucleotide microarray chip 
is a platform that was designed to include nearly all avail- 
able maize genes (Skibbe and Walbot, 2009). It includes 
expressed sequence tags involved in flowering, fruit forma- 
tion, root and shoot development, and germination, as well 
as developing and mature leaves (Skibbe and Walbot, 2009). 
A prior study utilizing this platform found that unstressed 
mature Z. mays leaves expressed approximately 26 000 tran- 
scripts under greenhouse growth conditions (Casati and 
Walbot, 2008). This indicates that approximately 59% of 
transcripts on the array are expressed in adult Z. mays leaves. 
A study on the diurnal rhythms of the maize leaf trancrip- 
tome, using a 105K Agilent maize array designed for all 
maize genomic and transcript sequences, found that approxi- 
mately 44 000 transcripts were detectable in the tested maize 
leaves, equating to 42% of the total transcripts (Hayes et al. , 
2010). Since the majority of transcripts associated with pho- 
tosynthesis are maximally expressed during daylight hours 
(Urbanczyk-Wochniak et al, 2005), samples were taken in 
this study 2-3 h into the photoperiod. At this time point, 
48% of the probes were expressed in M. x giganteus leaves in 
the control (25 °C) and chilling (14 °C). This is comparable 
to other studies using Agilent maize arrays. This hybridiza- 
tion efficiency for M. x giganteus leaf mRNA to the maize 
array is consistent with the >90% sequence identity observed 
among exon sequences between M. x giganteus and Z. mays 
(Swaminathan et al, 2010). This is not surprising given that 
both are members of the same grass tribe, Andropogoneae. 
This grass tribe is exclusively of the NADP-ME C 4 subtype 
(Giussani et al, 2001), suggesting that C 4 photosynthesis 
in maize, Miscanthus, and other species of the tribe most 
likely evolved from a single common ancestor (Sage, 2004; 
Swaminathan et al, 2010). 

While Z. mays is a crop of tropical origin, more than 60% 
of global production today is in the temperate zone (Long 
and Spence, 2013). It is therefore grown in colder climates 
more than any other C 4 food crop, yet its inability to accli- 
mate photosynthesis to chilling conditions is an Achilles heel 
(Baker et al., 1989; Dohleman et al, 2009). It is assumed 
that C 4 photosynthesis evolved in hot climates and subse- 
quently radiated to more temperate climates. The photosyn- 
thetic genes identified as upregulated in chilling-tolerant M. 
x giganteus in this comparative study indicate targets for 
upregulation, by selection or bioengineering, to improve the 



chilling tolerance of maize and to contribute to allowing a 
longer growing season for the crop in temperate climates. 

Supplementary material 

Supplementary data are available at JXB online. 

Supplementary Table SI. All transcripts of M. x giganteus 
leaves found to be significantly upregulated by treatment with 
a log 2 fold-change >0.70. 

Supplementary Table S2. All transcripts of M. x giganteus 
leaves found to be significantly downregulated by treatment 
with a log 2 fold-change >0.70. 

Supplementary Table S3. Transcripts for Rubisco, Rubisco- 
interacting proteins, and PPDK transcripts showing no significant 
log 2 fold-change and maintaining high expression during chilling. 
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